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BACKGROUND OF THE INVENTION 

15 Field of the Invention 

The present invention relates to a heat treatment apparatus , 
a controller for the heat treatment apparatus and a method 
for controlling the heat treatment apparatus, more 
specifically to a heat treatment apparatus in which 

20 temperatures of objects-to-be-processed can be 
incontiguously estimated, and a controller for the heat 
treatment apparatus and a method for controlling the heat 
treatment apparatus . 

25 R e l ated Backg rou nd Art 

One apparatus for making heat treatments on objects- 
to-be-processed, such as semiconductor wafers (hereinafter 
called wafers) or others, which is used in semiconductor 
device fabrication methods is a vertical heat treatment 

30 apparatus for batch treatment. in this vertical heat 
treatment apparatus, a number of wafers are held horizontally 
in a shelves-like manner on a holder, such as wafer boat or 
others, and the holder is loaded in the vertical heat 
treatment apparatus for a heat treatment, e.g. , CVD (Chemical 

35 Vapor Deposition), oxidation, etc. 

To make a heat treatment on the wafers it is necessary to 
precisely control temperatures of the wafers. For example, 
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in forming a thin film on the wafers by, e.g., CVD, a film 
thickness of the thin film depends on temperatures of the 
wafers. Accordingly, the temperature control of the heat 
treatment apparatus must be performed with high accuracy. 

5 The heat control has been conventionally made by loading 
the wafers with thermocouples into the heat treatment 
apparatus to meter temperatures of the wafers. 

However, when the wafers with thermocouples are loaded in 
the heat treatment apparatus, there is a risk that the metal 

10 forming the thermocouples disperses in the heat treatment 
apparatus and stays on the inside of the heat treatment 
apparatus, with a result that the staying metal adheres to 
the wafers, causing metal contamination. 

15 STTMMATRV Of INVENTION 

The present invention was made in view of the above 
circumstances, and an object of the present invention is to 
provide a heat treatment apparatus , a controller for the heat 
treatment apparatus and a method for controlling the heat 

20 treatment apparatus which are free from the risk of 
contamination of objects-to-be-processed and can control 
temperatures with high accuracy. 

(1) To attain the above-described object, the heat 
treatment apparatus according to the present invention 

25 comprises: a heater for heating an object-to-be-processed; 
a temperature meter disposed in the heater; a temperature 
estimator for computing an estimated temperature of the 
object-to-be-processed, based on a temperature metered 
result of the temperature meter; an error estimator for 

30 computing an estimation error of the estimated temperature 
computed by the temperature estimator; a temperature 
corrector for computing a corrected estimated temperature 
given by correcting the estimated temperature computed by 
the temperature estimator, based on the estimation error 

35 computed by the error estimator; and a heating controller 
for controlling the heater, based on the corrected estimated 
temperature computed by the temperature corrector, and a 



temperature recipe stating a relationship between a set 
temperature and a time. 

The temperature meter disposed in the heater meters 
temperatures, and, based on a temperature metered result, 
a temperature of the object-to-be-processed are estimated. 
The estimated temperatures are corrected, based on estimation 
errors. Consequently, a temperature of the object-to-be- 
processed can be accurately estimated. 

(2) The controller for controlling a heat treatment 
apparatus according to the present invention comprising a 
heater for heating an object-to-be-processed, and 
temperature meter disposed in the heater comprises a 
temperature estimator for computing an estimated temperature 
of the object-to-be-processed, based on a temperature 
metering result of the temperature meter; an error estimator 
for computing an estimation error of the estimated 
temperature computed by the temperature estimator; a 
temperature corrector for computing a corrected estimated 
temperature given by correcting the estimated temperature 
computed by the temperature estimator, based on the 
estimation error computed by the error estimator; and a 
heating controller for controlling the heater, based on the 
corrected estimated temperature computed by the temperature 
corrector. 

(3) The controller for controlling a heat treatment 
apparatus comprising a heater for heating an object-to- 
be-processed comprises a heating controller for controlling 
the heater in accordance with a temperature recipe stating 
a relationship between a set temperature and a time; and a 
control period determining unit for determining a control 
period for periodically controlling the heater by the heating 
controller, based on a change rate of the set temperature. 

Periods of the control are changed, based on change rates 
of a set temperature, whereby both accuracy and efficiency 
of the control can be considered. 
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(4 ) The method for controlling a heat treatment apparatus 
comprising a heater for heating an object-to-be-processed 
comprises a temperature metering step of metering a 
5 temperature in the heater; a temperature estimating step of 
computing an estimated temperature of the object-to-be- 
processed, based on a temperature metering result of the 
temperature metering step; an error estimating step of 
computing an estimation error of the estimated temperature 
10 computed in the temperature estimating step; a temperature 
correcting step of computing a corrected estimated 

0 temperature given by correcting the estimated temperature 
Jl computed in the temperature estimating step, based on the 
9^ estimation error computed in the error estimation step; and 

2 m a heating control step for controlling the heater, based on 

03 the corrected estimated temperature computed in the 

1 temperature correcting step, and a temperature recipe stating 
O a relation ship between a set temperature and a time. 

m 

SJ 20 The temperature meter disposed in the heater meters 

temperatures, and, based on temperature metered results, 
temperatures of the object-to-be-processed are estimated, 
and the estimated temperatures are corrected based on 
estimation errors. The heater is controlled, based on the 
25 corrected estimated temperatures . Consequently, the heater 
can be controlled based on the estimated temperatures of the 
object-to-be-processed computed with good accuracy. 

( 5 ) The method for controlling a heat treatment apparatus 
30 comprises a heating control step of controlling the heater 
in accordance with a temperature recipe stating a 
relationship between a set temperature and a time; and a 
control period determining step of determining a control 
period for periodically controlling the heater by the heating 
35 control step, based on a change rate of the set temperature. 

Periods of the control are changed, based on change rates 
of a set temperature, whereby both accuracy and efficiency 



of the control can be considered. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a partial sectional view of the heat treatment 
5 apparatus according to the present invention, which shows 
a structure thereof. 

FIG. 2 is perspective views of the heat treatment apparatus 
according to the present invention, which show the structure 
thereof. 

10 FIG. 3 is a block diagram of the controller of the heat 

treatment apparatus according to the present invention. 

FIG. 4 is a flow chart of a control procedure of the 
controller. 

FIG. 5 is a table which is one example of the temperature 
15 recipe. 

FIG. 6 is a graph which is another example of the 
temperature recipe. 

DES CR IPT ION O F THE PREFERRED EM BOD IMENT S 
20 (A First Embodiment) 

The vertical heat treatment apparatus according to a first 
embodiment of the present invention will be explained. 

FlGs. 1 and 2 are respectively a partial sectional view 
and perspective views of the vertical heat treatment 
25 apparatus according to the present embodiment. 

The vertical heat treatment apparatus according to the 
present embodiment comprises a reaction tube 2 of the 
double-tube structure of an inner tube 2a and an outer tube 
2b which are formed of, e.g., quartz, and a cylindrical 
30 manifold 21 of metal disposed on the bottom of the reaction 
tube 2. 

The inner tube 2a has the top uninterrupted and is supported 
on the inside of the manifold 21. The outer tube 2b has the 
top interrupted and has the lower end connected air-tightly 
35 to the upper end of the manifold 21 below a base plate 22. 

As shown in FIG. 2, in the reaction tube 2, a number of 
objects-to-be-processed, semiconductor wafers W are held on 



a holder, a wafer boat 23 horizontally and vertically spaced 
from each other in a shelves-like manner. The wafer boat 23 
is retained on a lid 24 through a heat insulation cylinder 
(heat insulator) 25. 
5 A heater 3 for heating the wafers W, objects-to-be- 
processed is disposed around the reaction tube 2, and the 
heater 3 is in the form of, e.g., resistance heaters. The 
heater 3 is divided in zone 1 to zone 5. The respective 
heaters 31-35 can have calories controlled independent of 

10 each other. Electric power controllers 41 - 45 have 
respective electric power sensors Sip - S5p for metering 
thermal powers of the respective heaters 31-35. 

Inside temperature sensors Slin - S5in in the form of 
thermocouples or others as an objects-to-be-processed 

15 vicinity temperature meters are disposed on the inside wall 
of the inner tube 2a corresponding to the zones 1 - 5 of the 
respective heaters 1 - 5 . Outside temperature sensors Slout 

- SSout in the form of thermocouples or others as a heater 
vicinity temperature meters are disposed on the outside wall 

20 of the outer tube 2b corresponding to the zones 1 - 5 of the 

respective heaters 31 - 35. 

Monitor wafers Wl - W5 are mounted at positions 

corresponding to the zones 1 - 5 of the respective heaters 

31-35. As will be described later, temperatures of the 
25 monitor wafers Wl - W5 are estimated, based on metered signals 

of the temperature sensors Sin (Slin - S5in), Sout (Slout 

- SSout) and the electric power sensors Sp(Slp -S5p) . 
The manifold 21 has a plurality of gas feed pipes for 

feeding gases into the inner tube 2a. In FIG. 1, two gas feed 
30 pipes 51 , 52 are shown for the convenience of the description. 
Flow rate controllers 61, 62, such as mass flow controllers 
or others, for controlling respective gas flow rates, and 
valves (not shown) are inserted in the gas feed pipes 51, 
52. 

35 The manifold 21 is connected to an exhaust pipe 27 for 

exhausting the gases through a gap between the inner tube 
2a and the outer tube 2b. The exhaust pipe 27 is connected 
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to a vacuum pump. A pressure adjuster 28 including a 
butterfly valve and a valve operator, for adjusting pressures 
in the reaction tube 2 is inserted in the exhaust pipe 27. 
The vertical heat treatment apparatus comprises a 
5 controller 100 for controlling treatment parameters, such 
as temperatures of a treatment atmosphere in the reaction 
tube 2, pressures in the reaction tube 2 and gas flow rates. 
The controller 100 receives metered signals of the electric 
power sensors Sp (Sip - S5p) and the temperature sensors Sin 
10 (Slin - S5in), Sout (Slout - S5out) and outputs control 
signals to the electric power controllers 41 - 45 of the heater 
3, a pressure adjuster 28 and flow rate adjusters 61 f 62. 
Then, the controller 100 will be described in good detail. 
FIG. 3 is a block diagram of a part of an internal structure 
15 of the controller 100, which is involved in the control of 
the heater 3. 

As shown in FIG. 3, the controller 100 comprises an A/D 
(Analog/Digital) converter 110 which converts analog metered 
signals of the electric power sensors Sp and the temperature 
20 sensors Sout, Sin to digital metered signal; a temperature 
estimator 120 which computes estimated temperatures of 
wafers; an error estimator 130 which computes estimation 

errors AT of the wafer estimated temperatures T' computed 
by the temperature estimator 120; an adder 140 as a 
25 temperature corrector which adds the computed estimated 

temperatures T' and the estimation errors AT to each other 
to compute corrected estimated temperatures T' ' of the 
wafers; a temperature recipe storage 150 which stores a 
temperature recipe stating relationships between set 

30 temperatures and times for a heat treatment of the wafers; 
a heating controller 160 which outputs electric power control 
signals P' to the electric power controllers 41-45, based 
on the corrected estimated temperatures T 7/ and the set 
temperatures Tsp set on the temperature recipe stored in the 

35 temperature recipe storage 150; and a computation period 
determining unit which determines a computation period ts 
for the computation of the temperature estimator 12 0 and the 
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error estimator 130. 

The temperature estimator 120 includes a first data 
interpolator 122 for interpolating data necessary to compute 
estimated temperatures T' of the wafers W corresponding to 
5 a computation period ts, and a temperature estimation model 
module 124 which computes estimated temperatures T ' of the 
wafers. The error estimator 13 0 includes a second data 
interpolator 132 which interpolates data necessary to compute 
estimation errors AT of the estimated temperatures of the 
10 wafers W, and an error estimation model module 134 which 
computes estimation errors AT of the estimated temperatures 
of the wafers. 

FIG. 4 is a flow chart of a control procedure of controlling 
15 the heater 3 . The procedure of the temperature control will 
be explained with reference to FIG. 4. 

(A) The computation period determining unit 17 0 determines 
a computation period ts, based on a change ratio RR of a set 
20 temperature Tsp stated in a temperature recipe stored in the 
temperature recipe storage 150 (Sll). 

A computation period ts is set corresponding to a change 
ratio RR of a set temperature, i.e. , a change ratio of a wafer 

25 temperature, whereby the computing and the control devices 
can be effectively used. That is, when a change ratio RR of 
a set temperature Tsp is small, a computation interval can 
be large, which allows the computing or control elements, 
such as CPU (Central Processing Unit) not shown to be used 

30 for other processing. On the other hand, when a change ratio 
RR of a set temperature is large, estimated temperatures, 
etc. of the wafers are computed by a short computation period, 
and temperatures of the wafers can be accordingly precisely 
controlled. 

35 

The temperature recipe stored in the temperature recipe 
storage 150 shows relationships between set temperatures Tsp, 
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i.e. f target temperatures of a heat treatment for the wafers 
W f and time. FIGs. 5 and 6 show examples of the temperature 
recipe. 

5 FIG. 5 shows a table expressing set temperatures Tsp (Set 
Point) and change rate RR of the set temperatures Tsp 
corresponding to sections of a time. FIG. 6 is a graph 
expressing the set temperatures Tsp and the sections of a 
time. 

10 

FIGs. 5 and 6 are different from each other in the forms 
and express the same temperature recipe. From a time to to 
a time tl, a set temperature is kept constant at 300°c, and 
from a time tl to a time t2, the set temperature Tsp changes 

15 from 300°C to 310°C at a 10 °C/min change rate. From a time 
t2 to a time t3, the set temperature Tsp changes respectively 
from 310°C to 350°C / from 350°C to 400°C and from 400° to 500°C 
respectively at a 20°C/min change rate RR, a 50°C/min change 
rate RR and a 100°C/min change rate RR. From a time t5 to 

20 a time t6, the set temperature is kept constant at 500°C. From 
a time t6 to a time t7 , from a time t7 to a time t8 and from 
a time t8 to a time t9, the set temperatures lower respectively 
from 500°C to 400°C, from 400°C to 350°C and from 350°C to 300°C 
respectively at a -100°C/min change rate, a -50°C/min change 

25 rate and a -10°C/min change rate. From a time t9 later, the 
set temperature is kept constant at 300°C. 

FIGs. 5 and 6 are substantially the same, and either of 
the forms may be used in practicing the present invention. 
30 In other words, the temperature recipe may be expressed in 
any form which can give set temperatures Tsp and change rates 
RR of the set temperatures corresponding to times. 

Based on change rates RR of set temperatures, computation 
35 periods ts are determined. One example of the method for 
determining set temperatures is expressed by the following 
formulae ( 1 ) to ( 4 ) . 
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ts=0.5[sec] (|RR|^50[°C/min] ) ... (1) 

ts=1.0[sec] (25s|RR|<50[°C/min] ) ... (2) 

ts=2.0[sec] (20£|RR|<25[°C/min] ) ... (3) 

5 ts=4.0[sec]( 0<s|RR|<10[°C/min] ) ... (4) 

Formula (1) gives a shortest computation period At=0.5 sec 
when a change rate RR of a set temperature is above 50°C/min 
including 50°C/min. The computation period ts changes from 
10 0.5 sec to 4.0 see, becoming smaller as a change rate of an 
absolute value of a set temperature becomes larger. 

In the example shown by Formulae (1) to (4), computation 
periods ts have values given by multiplying a shortest 
computation period by powers of 2. However, a computation 
period ts may be integer times a shortest computation period 
At. 

(B) The data interpolator 122 judges whether or not a 
20 computation period ts is equal to a shortest computation 
period At (S12) . 

When the judgement is YES, the temperature estimation model 
module 124 computes an estimated temperatures T' of the wafers 
25 (S13 ) . When the judgement is NO in S12 , the data interpolator 
122 interpolates data necessary for the computation of the 
temperature estimation (S14). The temperature estimation 
model module 124 computes an estimated temperature T' (S13) . 

30 To compute an estimated temperatures T of the wafers, 
digital signal outputs of the A/D converter 110, which include 
wafer vicinity temperatures Tin (Tlln - T5in), heater 
vicinity temperatures Tout (Tlout - TSout) and heat powers 
of the haters (PI - P5) (control signals of the controller 

35 100 for the electric power controllers 41 - 45), are used. 
The A/D converter 110 periodically converts analog signal 
inputs to digital signal outputs to output the digital signals . 
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To effectively use the computing and control devices 
including the A/D converter 110 it is preferable that a 
sampling interval agrees with a shortest computation period 

At. 

1) Here, the temperature estimation in S13 will be 
detailed 

Here, estimated wafer temperatures T' (estimated central 
temperatures TC (Tel' - Tc5') in a vicinity of the centers 
of the wafers) , estimated edge temperatures Te' in a vicinity 
of the edges of the wafers (Tel' - Te5 ' ) and estimated wafer 
vicinity temperatures Tin' (Tinl' - Tin5') are calculated, 
based on wafer vicinity temperatures Tin, heater vicinity 
temperatures Tout and heat powers of the haters (Pi - P5) 
(control signals of the controller 100 for the electric power 
controllers 41 - 45). The subscripts 1 to 5 correspond 
respectively to zones 1 to 5. 

y(k)=-al-y(k-l)-a2-y(k-2)-a3-y(k-3)- • -ai-y(k-i)-- • - 
al6-y(k-16)+bl-u(k-l)+b2-u(k-2)+b3-u(k-3)+- • -+bi-u(k- 
i) + - • -+bl6-u(k-16)+w(k) ... (5) 

wherein y(k)is an output vector after one period (At [sec] 
later); y(k-l) is a current output vector; y(k-i) is an 
output vector before (i-1) period (before (i-1) -At[sec] ) ; 
u(k-l) is a current input vector; u(k-i) is an input vector 
before (i-1) period (before (i-1) -At[sec] ) , and w(k) is a 
noise vector (white noise) . 

vectors u(k-i) , y(k-i) andw(k) are specifically expressed 
by the following formulae ( 6 ) to ( 8 ) . 

u(k-i) = (PI, P2, P5, 

Tlout, T2out, T6out, 

Tlin, T2in, T5in) ... (6) 

y(k-i) = (TIC, T2C, T5c ' , 
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12 

Tle f f T2e', T5e', 

Tlin', T2in', T5in' ) .« (7) 

w(k) = (wl(k), w2(k), wl5(k)) »• (8) 

5 

Formula (5) is a 16 dimension-parametric model, and 
computes an output vector y(k) after one period, based on 
a current and a past input vectors u(k-i) and a current and 
a past output vectors y(k-i). 

10 

An output vector y(k) after one period is used as a current 
value y(k-l) of an output vector upon the next computation 
(at which 1 is added to k) . When an input vector u(k-l) is 
known, an output vector y(k) after one period can be 
15 sequentially computed. 

That is, estimated wafer central temperatures Tc'(Tlc f - 
T5C), estimated wafer edge temperatures Te' (Tie' - T5e') 
and estimated wafer vicinity temperatures Tin' (Tlin' - 
20 T5in') can be computed based on heat powers P (Pi - P5) of 
the haters, heater vicinity temperatures Tout (Tlout - T5out) 
and wafer vicinity temperatures Tin (Tlin - T5in) . 

The estimated wafer vicinity temperatures Tin' are 
25 computed for the error estimation in the later step Si 6 • When 
the error estimation is unnecessary, the output vector y (k-i) 
may not include an estimated wafer vicinity temperatures Tin' 
as in the following Formula (9). 

30 y(k-i) = (Tlc' , T2c ' , T5C 

Tie, T2e, T5e) — (9) 

In this case, both estimated wafer temperatures T' and 
estimated wafer vicinity temperatures Tin' are computed by 
35 one model but may be computed respectively by respective 
models. For example, estimated wafer temperatures T' are 
computed by a model using the following formulae ( 10 ) to ( 13 ) , 
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and estimated wafer vicinity temperatures Tin' are computed 
by a model using the following formulae (14) to (17). 

yl(k)=-all-yl(k-l)-al2-yl(k-2) 

-• • --aliyl(k-i) -• • • -all6 -yl (k-16 ) 

+bllul(k-l)+bl2-ul(k-2) 
+• • -+bli-u(k-i)+ • • +bll6 ul(k-16) 

+wl(k) ••• (10) 

ul(k-i) = (Pi, P5, 

Tlout, ~, T5out, 

Tlin, T5in) ... (11) 

yl(k-i) = (TIC, T5c ' , 

Tie', T5e) .» (12) 

wl(k) = (wll(k), wll6(k)) - (13) 

y2(k) = -a21y2(k-l)-a22*y2(k-2) 

• --a2i-y2(k-i)-- • --a216-y2(k-16) 
+b21u2(k-l)+b22u2(k-2) 
+• • -+b2i-u2(k-i)+- • -+b216-u2(k-16)+w2(k) 

... ( 14 ) 

u2(k-i) = (PI, P5, 

Tlout, T5out) - (15) 

y2(k-i) = Tlin', T5in') - (16) 

w2(k) = (w21(k), w216(k)) ... (17) 

Parameters al - al6, bl - bl6 of Formula (5) and etc. and 
a noise vector w( k) are determined in advance before an output 
vector y(k) is computed. 



To determine parameters al - al6, bl - bl6, a noise vector 
w(k), subspace method or Auto-Regressive Exogeneous Model 
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(hereinafter called ARX model) can be used. 

Specifically, data, such as metered signals of the 
temperature sensors Slin - S5in, slout - S5out and real 
5 metered temperatures of the monitor wafers Wl - W5 
(temperature sensors, such as thermocouples are mounted on 
the monitor wafers) are inputted to, e.g. software Matlab 
(produced by The Math Works , Inc. , marketed by Cybernet System 
Kabushiki Kaisha) to inversely compute the parameters al - 
10 al6, bl - bl6, and the noise vector w(k). 

A plurality of combinations of the given parameters al - 
al6, bl -bl6andthe noise vector w(k) are usually available. 
Out of the plural combinations, one in which the estimated 
15 temperatures Tl ' - T5 ' given by Formula ( 5) well agree with 
the actually metered temperatures of the monitor wafers is 
selected (Model appreciation) . 

Dimensions of the model are usually large, and dimensions 
20 are suitably reduced to make the model of 16 dimensions. 

2) Next, the data interpolating step (S14) will be 
explained. 

25 When a computation period ts is equal to a shortest 

computation period At, a value of the output vector y(k- 
1) computed before can be substituted as it is into Formula 
(5) . However, when a computation period ts is different from 

(larger than) the shortest computation period At, the data 
30 is insufficient to be substituted as it is into Formula (5) 
or others. Then, the data must be interpolated. 

Here, it is assumed that a computation period ts is n- 

times a shortest computation period :t (ts=n*At). In this 
35 case, the computation of an output vector y(k) by Formula 
(5) is performed every n-shortest computation periods. That 
is, only when k=n-m in Formula (5), the computation is 



15 

performed (m: an integer). 

When n=4, for example, k is given by Formula (18). The 
data of the output vector with, e,g., k=l to 3 and 5 to 7 
5 must be interpolated. 

k=0, 4, 8, 12, 16, _ (18) 

In Formula (5) , when k=n*m, an output vector y(k) computed 
10 immediately before is k=n*(m-l), and an output vector 
y(k)computed next immediately before is k=n* (m-2). In such 
case, the following interpolation is performed. 

(1) Data interpolation between k=n-m and k=n*(m«l) 
15 First, an interpolation for an output vector between k=n-m 

and k=n* (m-1) will be explained. Output vectors y(k) for 
k=(n*m-l), (nm-2), (n* (m-l)+l) must be given. All the 
output vectors can have a value of the immediately before 
computed output vector (n*(m-l)) as expressed by the 
20 following Formula (19). 

y(n*m-l) = y(n* (m-1)) 
y(n-m-2) = y(n- (m-1)) 



y(n*(m-l)+l) = y(n-(m-l)) ... (19) 

25 

(2) Data interpolation between k-n-(m-l) and k=n*(m- 
2) output vectors 

An output vector between k=n- (m-1) and k=n* (m-2) can be 
an interpolated value between an output vector y(n*(m-l)) 

30 computed immediately before and an output vector y(n* (m- 
2) computed next immediately before. As exemplified by the 
following Formula (20), the interpolation can use linear 
approximation between the immediately before output vector 
y(n- (m-1)) and the next immediately before output vector 

35 y(n- (m-2)). 



y(n- (m-l)-l) 
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- [l-y(n- (m-2)) + (n-1) -y(n- (m-1 ) ) ]/n 
y(n* (m-l)-i) 

= [i*y(n- (m-2) ) + (n-i) -y (n- (m-1 ) ) ] /n 

5 

y(n* (m-2)+l) 

= [(n-l) -y(n-(m-2)) + l-y(n*(m-l)) ]/n ... (20) 

As the interpolation, various interpolations, such as 
10 parabolic approximation beside the linear approximation 
expressed by Formula (20) can be used. 

The data interpolation between k=n-m and k=n- (m-1) can be 
made by interpolations besides the interpolation expressed 
15 by Formula (19) (a computed value of an output vector y(k) 
computed immediately before). For example, extrapolation 
maybe used, based on a value of an output vector computed 
immediately before and a value of an output vector next 
immediately before. 

20 

The following Formula (21) exemplifies linear 
extrapolation using values of output vectors computed 
immediately before and next immediately before. 

25 y(n*(m-l)) 

= [<2n-l) -y(n'(m-l)) 

-(n-1) -y(n*(m-2))]/n 



y(n- (m-l)+i) 
30 = [(n+i) -y(n- (m-1) )-i-y(n* (m-2) ) ]/n 



y(n- (m-l)+l 

= [(n+1) •y(n-(m-l))-l-y(n-(m-2))]/n ~ (21) 



35 



(3) Data interpolation before k=n*(m-2) 
A value of an output vector y(k-i) before k=n*(m-2) may 
be value of the value of the output vector y(k-i) used in 
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computing the output vector y(k) immediately before. 

Data are thus interpolated by using Formulae (20) and (21) 
for the temperature estimation. 

5 

In the above-described interpolation, when a computation 
period is switched from a shorter computation period to a 
longer computation period, an output vector y(k-i) computed 
in the shorter computation period can be used in the 

10 computation of Formula (5) early after the switch of the 
computation period- In short, when a suitable value of the 
output vector y(k-i) is available, the value is used, and 
when a suitable value is unavailable, the interpolation is 
performed by an interpolation or an extrapolation to give 

15 a current or a past output vector. 

(C) It is judged whether or not a computed period is equal 
to a shortest computation period ( S15 ) , and when the judgement 
is YES, estimation errors AT of wafer temperatures are 
20 computed (S16). When the judgement in S15 is NO, data 
necessary for the computation of the temperature estimation 
are interpolated (S17), and then estimation errors AT are 
computed (16). 

25 The error estimation in S15 will be detailed. 

The error estimation is for estimating differences between 
values of estimated wafer temperatures T' (Tc f , Te') and 
actually metered temperatures T (Tc, Te). That is, 
30 idealistically, estimation errors AT (central estimation 
errors ATc in a vicinity of the centers of the wafers and 
edge estimation errors ATe) are expressed by the following 
Formula (22). 

35 AT = T - T' (22) 



Actual temperatures T of the wafers are not actually 
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metered, and estimation errors AT of the wafers must be 
estimated, based on relationships with other parameters. As 
such parameter, estimated wafer vicinity temperatures Tin f 
and actually metered wafer vicinity temperatures Tin can be 
used. Here, an estimated wafer vicinity estimation error 
ATin can be given by the following Formula (23). 

ATin = Tin - Tin' ... (23) 

It is considered that temperatures of the wafers themselves 
and temperatures in the vicinity of the wafers have close 
relationships with each other. Empirically, estimation 
errors AT have values related with vicinity estimation errors 
ATin. Accordingly, a model is prepared to thereby compute 
estimation errors AT of the wafers, based on vicinity 
estimation errors ATin. Change rates RR of set temperatures 
have close relationship with estimation errors AT of the 
wafers, and are taken into consideration. 

Finally, estimation errors AT of the wafers are computed, 
based on vicinity estimation errors ATin and change rates 
RR of set temperatures. 

To compute estimation errors in S16, a parameter model 
expressed by the following Formula (24) which is 
substantially the same as Formula (5) can be used. 

y(k) = -al-y(k-l)-a2-y(k-2)-a3*y(k-3)- • * • -ai-y (k-i) 
- *-al6-y(k-16) + bl * u ( k-1 )+b2 • u ( k-2 ) 
+ b3*u(k-3)+ • • -+bi*u(k-i) + *♦ * 
+bl6-u(k-16) + Aw(k) (24) 

wherein y(k) is an output vector after one period (At [sec] 
later); y(k-l) is a current output vector ; y(k-i) is an output 

vector before (i-1) period ((i-1) -At [sec] before); u(k- 
1) is a current input vector; u(k-i) is an input vector before 
(i-1) period ((i-1) *At[sec]); and w(k) is a noise vector 
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(white noise), which are substantially the same as those in 
Formula ( 5 ) . 

However , vectors u(k-i), y(k-i) andw(k) are specifically 
5 expressed by the following Forraulae(25) to (27), which are 
different from those of Formula (5). 

u(k-i) = (ATlin, AT2in, AT5in, RR) ... (25) 

10 y(k-i) = (ATlc, AT2c, AT5c, ATle, AT2e, AT5e) 

... (26) 

Aw(k) = (Awl(k), Awl5(k)) „. (27) 

15 Formula (24) is a 16 dimension-parametric model, as is 

Formula (5) . When an input vector is known, an output vector 
y(k) after one period can be sequentially computed* For this 
computation, parameters aO - al5, bO - bl5, and a noise vector 

Aw(k) can be determined by subspace method or ARX model, as 
20 can be used in Formula ( 5 ) . 

The data interpolating step (S17) can be carried out by 
an interpolation and an extrapolation substantially in the 
same way as in S14. The interpolation and the extrapolation 
25 are not detailed not to repeat their explanation. 

D) Temperature correction is performed, based on a 
temperature estimation result of S13 and a result of the error 
estimation in S16 (S18). 



30 



35 



Corrected estimated temperatures T ' ' can be given by adding 

estimation errors AT to estimated temperatures T' of the 
wafers by the adder 140 as the corrector as expressed by 
Formulae (28) and (29) • 

Tic" = Tic' + ATic (i = 1 - 5) ... (28) 

Tie" = Tie' + ATie (i = 1 - 5) ... (29) 



20 

Resultantly, temperatures of the wafers can be more 
correctly estimated* 

5 (E) The heater is controlled based on corrected estimated 

temperatures T' f (S19). 

The heater controller 160 compares corrected estimated 
temperatures T' ' with a set temperature Tsp to determined 
10 suitable heat powers PI' - P5' to the heaters 31 - 35. 

For example, powers of the heaters 31-35 can be determined 
corresponding to differences between corrected estimated 
temperatures T ' ' and a set temperature Tsp. This control is 
15 made corresponding to a computation period ts determined by 
the computation period determining unit 170. 

(Other Embodiments) 

The above-described embodiments of the present invention 
20 can be expanded and modified in the scope of the technical 
idea of the present invention ♦ 

The heat treatment apparatus according to the present 
invention does not essentially include both of the 
computation period determining unit and the temperature 
25 estimator, and may include either of them. For example, the 
computation period determining unit can be used to actually 
meter wafer temperatures, and in this case, the temperature 
control is effectively performed, based on computation 
periods determined by the computation period determining 
30 unit. 

The model is not essentially the parameter model expressed 
by Formula (5) and may be other suitable models. Dimensions 
of the model are not essentially 16 and may be 8, 20 or others. 

35 

The heater may not be divided and is not essentially divided 
in 5 sections. 



21 



The heat treatment apparatus according to the present 
invention is not essentially a vertical heat treatment 
apparatus or of the batch-type, and may be of sheet-type for 
making a heat treatment on wafers sheet by sheet. 

As described above, according the present invention, 
temperatures of objects-to-be-processed can be estimated 
with good accuracy. 



